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Elementary excitation in a supersolid 
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We study elementary low energy excitations inside a supersolid. We find that the coupling between 
the longitudinal lattice vibration mode and the superfluid mode leads to two longitudinal modes ( 
one upper branch and one lower branch ) inside the supersolid, while the transverse modes in the 
supersolid stay the same as those inside a normal solid. We also work out various experimental 
signatures of these novel elementary excitations by evaluating the Debye-Waller factor, density- 
density correlation, vortex loop-vertex loop interactions, specific heat and excess entropy from the 
vacancies per mole. 



1. Introduction. A supersolid is a state which has 
both a solid order and a superfluid order. The possibility 
of a supersohd phase in "^He was theoretically speculated 
in 1970A. Over the last 35 years, a number of exper- 
iments have been designed to search for the supersolid 
state without success. However, recently, by using tor- 
sional oscillator measurement, a PSU group lead by Chan 
observed a marked 1 - 2% NCRI of solid ^He at - 
in bulk '^He^. The authors suggested that the supersolid 
state of '^He maybe responsible for the NCRI. The PSU 
experiments rekindled extensive both theoretical^i^i^i^'^ 
and experimenta l^' interests. The torsional oscilla- 
tor measurement is essentially a dynamic measurement, 
as suggested in^i, other possibilities can also lead to the 
NCRI observed in the Kim-Chan experiments. Obvi- 
ously, many other thermodynamic and equilibrium mea- 
surements are needed to make a definite conclusion. 
Some interesting physics near the finite temperature nor- 
mal solid to supersolid transition was explored in5.. The 
fundamental questions to be addressed in this paper is 
what are the elementary excitations in a supersolid at 
zero and very low temperature and what are the phys- 
ical consequences of these elementary excitations which 
can be probed by various experimental techniques such 
as X-ray scattering, neutron scattering, acoustic wave at- 
tenuation and heat capacity in solid Helium 4. In prin- 
ciple, if these elementary low energy excitations can be 
detected by these experiments can prove or disprove the 
existence of the supersolid in Helium 4. 



2. Elementary excitations in a SS. Classical non- 
equilibrium hydrodynamics in SS was investigated for a 
long timelii?-. These hydrodynamics break down at very 
low temperature where quantum fluctuations dominate. 
However, the quantum nature of the excitations in the 
SS has not been studied yet. Here, we will study the 
quantum characteristics of low energy excitations in the 
SS. For example, how the phonon spectra in the SS differ 
from that in a NS and how the SF mode in the SS differs 
from that in a SF. The following effective Lagrangian in 
the imaginary time describing the low energy excitations 



inside a SS can be seen just from symmetry point of view: 

+ ^[f^idrOf + plpdaOdfjO] + aaf)Uai3idr9 (1) 

where Ua = Ua{x,T) is the lattice displacement, the first 
term is the phonon part, the second term is the superfluid 
part, the last term is the crucial coupling between the 
phonon part and the Berry phase term from the super- 
fluid part which comes from aapUafsip^ drip term in the 
phase representation of the superfluid order parameter 
ip ~ e'^. pn is the normal density, Uaf3 — ^{daUp + dpUa) 
is the linearized strain tensor, Xa/s-yS are the bare elastic 
constants dictated by the symmetry of the lattice, it has 
5 (2) independent elastic constants for a hep ( isotropic 
) lattice. 6 is the phase and k is the SF compressibil- 
ity. For an hep crystal, p^p is the SF stiffness which has 
the same symmetry as aaf3 — azUanp + a±(Sai3 — Uanp) 
with n a unit vector points along the preferred axis of the 
crystal. For an isotropic or a cubic crystal a^^ = aS^f^. 
In the following, we discuss two extreme cases: isotropic 
solid and hep lattice separately. Usual samples are be- 
tween the two extremes, but can be made very close to 
an hep crystal. 

(a) Isotropic solid: A truly isotropic solid can only 
be realized in a highly poly-crystalline sample. Usual 
samples are not completely isotropic. However, we expect 
the simple physics brought about in an isotropic solid 
may also apply qualitatively to other samples which is 
very poly-crystalline. 

For an isotropic solid, Xa(3-fS — XSa/iS-^s + p{5a-y5ps + 
Sa&Sp-i) where A and v are Lame coefficients, ^ = 
P^Sa,f3, aa,i3 = a5a,j3- In {q,ijjn) spacc, Eqn[T] becomes: 

As = ]^[pni^l + {X + 2p)q^]\ui{q,uJn)\'^ 

-f aquJnUi{-q,-uJn)e{q,uJn) 
+ + /V]|tit(9,w„)p 



(2) 
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FIG. 1: The elementary low energy excitations inside a super- 
solid. The coupling between the phonon mode u)p = v-pq ( the 
upper dashed line ) and the superfluid mode uis = Vsq ( the 
lower dashed line ) leads to the two new longitudinal modes 
tj± = v±q { solid lines ) in the SS. The transverse mode stays 
the same as that in a normal solid and is not shown. 



where ui{q,ujn) = iqiUi{q,uJn) / q is the longitudinal com- 
ponent, Ut{q,ujn) = ieijqiUj(q,uJn)/q are transverse com- 
ponents of the displacement field. Note that Eqn[2]shows 
that only longitudinal component couples to the super- 
fluid 9 mode, while the two transverse components are 
unaffected by the superfluid mode. This is expected, be- 
cause the superfluid mode is a longitudinal density mode 
itself which does not couple to the transverse modes. 

From Eqn[21 we can identify the longitudinal- 
longitudinal phonon correlation function (uiui) and also 
(69) and {ui9) correlation functions. By doing the ana- 
lytical continuation iujn uj + i6, we can identify the two 
poles of all the correlation functions at uj^. = v'^q'^ where 
the explicit expressions of the two velocities v± are given 
in^, but are not needed in our discussions. It is easy to 
show that > Vp > Vs > and v'^ + v"^ > + v'^, 
but v+V- — VpVs, so + > Vp + Vs { see Fig.l ). 

If setting a — 0, then obviously, reduces to the 
longitudinal phonon velocity vf^ = (A -I- 2/i)/p„ and the 
superfluid velocity = ps/ k respectively. Of course, the 
transverse phonon velocity v^p = p/ Pn is untouched. For 
notation simplicity, in the following, we just use Vp for vip. 
Inside the SS, due to the very small superfluid density ps, 
it is expected that Vp > Vg- In fact, in isotropic solid He'^, 
it was measured that vip ~ 450 — 500m/ s,Vt ~ 230 ~ 
320m/s and Vg ~ 366m/s near the melting curved. The 
size of the coupling constant a was estimated to be ~ 0.1 
from the slope of the melting curve^ii^. So u+ ( w_ ) are 
about 10% above ( below ) Vp { Vg )■ 

The two longitudinal modes in the SS can be under- 
stood from an intuitive picture: inside the NS, it was 
argued inii that there must be a diffusion mode of va- 
cancies in the NS. Inside the SS, the vacancies condense 
and lead to the extra superfluid mode. So the diffusion 
mode in the NS is replaced by the SF 9 mode in the 
SS. Its coupling to the lattice phonon modes lead to the 
elementary excitations in a SS shown in Fig.l. 



(b) hep crystal: Usual single hep crystal samples may 
also contain dislocations, grain boundaries. Here we ig- 
nore these line and plane defects and assume that there 
are only vacancies whose condensation leads to the su- 
perfluid density wave inside the supersolid'''^. 

For a uni-axial crystal such as an hep lattice, the action 

is: 

jO-hcp = ^[pnidrUa]'^ + Kii{ul^ + U^y) + 2Ki2UxxUyy 

+ 2(i^n - Ki2)uly + Kiiiul, + ul,)] 

+ ^[nid^fff + pl{d,9f + pU{d.9f + {dy9f)] 

+ [a^dzUz + a^{dxUx + dyUy)]idr9 (3) 

If q is along z direction, namely q^ ^ O^qx ~ qy ~ 0, 
then Eqn[3]simplifles to: 

^hcp = ]^[Pn^^l+ K:i'iql]\u^{q^,ujn)\'^ 
- iazqzUJnUz{-q2,-uJn)9{qz,uJn) 

+ \[Pn^l+KMql/'^]\ut{qz,^n)\^ (4) 

where |ut(g^,a;„)p = |ua:(g^,w„)p -I- |uy(g^,w„)p stand 
for the two transverse modes with the velocity — 
Kii/Apn- The superfluid mode only couples to the longi- 
tudinal Uz mode, while the two transverse modes Ux,Uy 
are decoupled. Eqn[4] is identical to Eqn[2] after the re- 
placement Uz — > u/,X33 — > A -|- 2/i, a, so vf^^ = 
Kss/pn- It was found that vip ^ 540m/s,ft ^ 250m/s 
when qis along the z direction^'*. Fig.l follows after these 
replacements. 

Similarly, we can work out the action in the xy plane 
where qz = 0,qx ^ 0, qy ^ 0. Then Uz mode is decoupled, 
only Ux,Uy modes are coupled to the superfluid mode: 

+ 2Ki2UxxUyy + 2{Kii - Ki2)uly] 

+ lKdr9r + pl{d^9r] 
+ a±daUaidT-9 

+ l[PnidrUzf+K4i/A{da,Uzf] (5) 

where a, (3 = x,y. By comparing EqnO with Eqnl2l 
we can see that Ku ^ A + 2p,, K12 A, so vf^^ — 
Kii/pn and all the discussions in the isotropic case can 
be used here after the replacements. Fig.l follows af- 
ter these replacements. Namely, only the longitudinal 
component in the xy plane is coupled to the 9 mode, 
while the transverse mode in the xy plane with veloc- 
ity v^^y — {Kii — Ki2)/2pn is decoupled. Obviously the 
transverse mode along z direction Uz mode with the ve- 
locity vfz = K44/4pn is also decoupled. Note that the 
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two transverse modes have different velocities. It was 
found that vip ~ Abbm/ s,Vtz ^ '255m/ s,vt^y ~ 225m/ s 
when q is along the xy plane^^. 

Along any general direction q, strictly speaking, one 
can not even define longitudinal and transverse modes, 
so the general action Eqnl3] should be used. Despite the 
much involved 4x4 matrix diagonization in Ux,Uy,Uz, 0, 
we expect the qualitative physics is still described by 
Fig.l. 

Recent inelastic neutron scattering (INS) did not de- 
tect any changes in atom kinetic energy ^PnidrUo/f' 
in the temperature range T — 70 mK — 400 mK^^ 
and atomic momentum distribution function n{k) with 
no = within T = 80mK — bOOmK--. these facts ex- 
clude the existence of SS at T > 70mK. It was known 
that INS is a very powerful tool to measure the phonon 
spectra in a normal solid (NS). We expect that if SS 
indeed exists, the INS can also be used to detect the pre- 
dicted elementary low energy excitation spectra in the 
SS shown in Fig.l. Namely, the INS should be able to 
map out the dispersion of two longitudinal modes U!± in 
Fig.l and the two transverse modes when T < T55. The 
neutron scattering cross-section from the uj±l and the 
spectral weight distribution between the two ±L modes 
will be calculatedii. 

3. Debye- Waller factor in the X-ray scatter- 
ing from the SS: It is known that due to zero-point 
quantum motion in any NS at very low temperature, the 
X-ray scattering amplitude /(G) will be diminished by a 
Debye- Waller (DW) factor e^3<3 where Ua is the 
lattice phonon modes in EqnlT] In EqulU if the coupling 
between the u and 9 were absent, then the DW factor 
in the SS would be the same as that in the NS. By tak- 
ing the ratio Iss{G)/Ins{G) at a given reciprocal lattice 
vector G, then this DW factor drops out. However, due 
to this coupling, the (m^) in SS is different than that in 
NS, so the DW factor will not drop out in the ratio. In 
this section, we will calculate this ratio and see how to 
take care of this factor when comparing with the X-ray 
scattering data. 

The density order parameter at the reciprocal lat- 
tice vector G is pq{x,t) = e^G-u{x,T)^ then {pg{x,T)) = 
g-iG,G, ^ The Debye-Waller factor: 

I[G) = |(pg(f,r)}|2 = ^-G.G,(u^{x,r)u,{x^r)) (g) 

where the phonon-phonon correlation function is: 

{uiUj) = {uiui)qiqj + {utUt){5i;j - cjiqj) (7) 

where qiq^ = . 

Then substituting EqnlT] into Eqn[6] leads to: 

a(G) = IssiG)/lNs{G) - e-iG^(^«')' (8) 

where (Au^), ^ {uf{x,T))ss - {uf {x , t)) n s and the 
transverse mode drops out, because it stays the same 
in the SS and in the NS. 



From Eqnl2l it is easy to see that (Am^); < 0, namely, 
the longitudinal vibration amplitude in SS is smal/er that 
that in NS. The a(G)(T) = g-^'^'^'^"")' > 1. This 
is expected, because the SS state is the ground state 
at T < Tss, so the longitudinal vibration amplitude 
should be reduced compared to the corresponding NS 
with the same parameters /3„,A,/i. It is easy to show 
that {Au^)i{T = 0) < and {Au^)i{T) - {Au^)i{T = 
0) ^ > 0. Namely, the difference in the ratio will 
decrease as as the temperature increases. Of course, 
when T approaches T55 from below, the difference van- 
ishes, the a(G) will approach 1 from above, the SS turns 
into a NS. The density-density correlation function was 
studied in^. Unfortunately, very recent INS^^ did not 
detect the predicted anomaly in the Debye-Waller factor 
within T = lAOmK — IK . This fact indicates the absence 
of SS when T > UOmK. 

4. Vortex loops in supersolid In section 3, we 
studied the low energy excitations shown in the Fig.l 
by neglecting the topological vortex loop in the phase 
6. Here, we will study how the vortex loop interac- 
tion in the SS differ from that in the SF. For simplic- 
ity, in the following, we only focus on the isotropic case. 
The formulations can be generalized to the hep case 
straightforwardly. We can perform a duality transforma- 
tion on Eqn[l] to the vortex loop representation in terms 
of the 6 components anti-symmetric tensor gauge field 
a^i, — — a^^ and the 6 components anti-symmetric ten- 
sor vortex loop current j^j^ = -^e^uXcrdxdaO due to the 
topological phase winding \n 0. It is the most conve- 
nient to choose the Coulomb gauge dattap = to get 
rid of the longitudinal component, then the transverse 
component is at = ieaP'y<laO.[3-i/<l- It can be shown that 
flfp = 2|aQ/3p. Then EquHJis dual to: 

Cy = ^[pnUjl + {X + 2p + a'^/K)q'^]\ui{q,uJn)\'^ 

1 2 

+ o(9^/'« + ^n/Ps)l«tP ^ 9^|aOQp 

2 Ps 
- aq^ / KUi[-q,-uJn)at{q,iOn) 

+ i27rjg^aoa + i2T:f^paap (9) 

where the transverse phonon mode Ut was omitted, be- 
cause it stays the same as in the NS as shown in Eqn[2l 
It is easy to see that only at has the dynamics, while 
aoa is static. This is expected, because although a^^ has 
6 non-vanishing components, only the transverse com- 
ponent at has the dynamics which leads to the original 
gapless superfluid mode = v'lq'^ . EqnlQ] shows that 
the coupling is between the longitudinal phonon mode 
ui and the transverse gauge mode at- The vortex loop 
density is jg^ — ^eap-ydpd^B and the vortex current is 
■?a/3 = ^<^ai3'y[do,d'y]0- Integrating out the aoa, we get 
the vortex loop density-density interaction: 

TrPs J J dxdyj^^{x,T)j--^—^j^^{y,T) (10) 
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Namely, the vortex loop density- density interaction in SS 
stays as 1/r which is the same as that in SF. Therefore, 
a single vortex loop energy and the critical transition 
temperature T^dxy^ is solely determined by the superfluid 
density ps independent of any other parameters in Eqnl2l 
except that the vortex core of the vortex loop is much 
larger than that in a superfluid''. The critical behaviors 
of the vortex loops close to the 3d XY transition was 
studied ini^. 

Integrating out the at , we get the vortex loop current- 
current interaction: 

27r^jQ/3(-'L --'^n)Da0.-,5{q,^n)j'^s{?^^n) (H) 

where Dap.-^siq.i^n) = (Sa-ySps - Sp^S^s - ^^^a-y - 
^Sf}S + ^Sc.s + ^S^s)Dt{q,u;n) where is 
the at propagator. Defining ADt{q,uJn) — Df^{q,uJn) — 
Df^{q,ujn) as the difference between the at propagator 
in the SS and the SF. For simplicity, we just give the ex- 
pression for the equal time ADt{x — x\t — 0) — j^gjg,y} 
where c is a positive constant if assuming SS has the same 
parameters k, Ps as the SF. 

5. Specific heat in the SS It is well known that at 
low T, the specific heat in the NS is C^^ = C^p^ + Cl^p^ + 
Can where C,^^ = ^ksi^g^f is from the longitudi- 
nal phonon mode and C[^^ = 2 x ^ksij^)^ is from 
the two transverse phonon modes, while Cyan is from 
the vacancy contribution. C-uan was calculated ii)^ by 
assuming 3 different kinds of models for the vacancies. 
So far, there is no consistency between the calculated 
Cjjan and the experimentally measured one^ii^. The spe- 
cific heat in the SF Cf^ = TT^sCig^f)^ is due to the 
SF mode 9. From EqnO we can find the specific heat in 
the SS: 

where C*^ stands for the contributions from the trans- 
verse phonons which are the same as those in the NS. 

It was argued the critical regime of finite temper- 
ature NS to SS transition in Fig. 1 is much narrower than 
the that of SF to the NL transition, so there should be a 
jump in the specific heat dXT — Tss- From Eqn[T2l it is 
easy to see that the excess entropy due to the vacancies 
is: 

(13) 



where AS > is dominated by the lower branch v- < 
vip in Fig.l. Using the molar volume vq ^ 20cm^ / mole 
of solid ^He and Tss ^ lOOmX, we can estimate the 
AS* per mole is ~ lO^^i? where R is th gas constant. 
This estimate is 3 orders magnitude smaller than that 
iT_}\3^ where the SS state was taken simply as the boson 
condensation of non- interacting vacancies. Our estimate 
is indeed consistent with recent experiment on specific 
heat^. 



6. Conclusions: In this paper, starting from the 
quantum Ginsburg-Landau theory developed inii^, we 
worked out the elementary excitations inside a super- 
solid. We found that the elementary excitations have 
two longitudinal modes liJ± = v±q shown in Fig.l. The 
transverse modes in the SS stays the same as those in the 
NS. The uj± are estimated to be 10% higher ( lower ) than 
the sound speed in the normal solid and the superfluid 
respectively. Then we calculated the experimental signa- 
ture of the two modes. We found that the longitudinal 
vibration in the SS is smaller than that in the NS ( with 
the same corresponding solid parameters ), so the DW 
factor at a given reciprocal lattice vector is larger than 
that in the NS. The density-density correlation function 
in the SS is weaker than that in the NS. By going to 
the dual vortex loop representation, we found the vortex 
loop density-density interaction in a SS stays the same 
as that in the SF ( with the same corresponding super- 
fluid parameters ), so the vortex loop energy and the 
corresponding SS to NS transition temperature is solely 
determined by the superfluid density and independent of 
any other parameters. The vortex current-current inter- 
action in a SS is stronger than that in the SF. The specific 
heat in the SS is given by the sum from the transverse 
phonons and the two longitudinal phonons and still shows 
~ behavior. The excess entropy due to the vacancies 
was estimated to be 3 order of magnitude smaller than 
the previous idea bose gas estimation. Comparison with 
very recent neutron scattering measurements are made. 
No matter if the SS exists in '^He, the results achieved 
should be interesting in its own right and may have ap- 
plications in other systems such as the possible excitonic 
supersolid in electron-hole bilayer systems^. 
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